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Abstract 

We explore the possibility of simultaneous determination of neutrino mass hierarchy and the CP 
violating phase by using two identical detectors placed at different baseline distances. We focus 
on a possible experimental setup using neutrino beam from J-PARC facility in Japan with beam 
power of 4MW and megaton (Mton)-class water Cherenkov detectors, one placed in Kamioka and 
the other at somewhere in Korea. We demonstrate, under reasonable assumptions of systematic 
uncertainties, that the two-detector complex with each fiducial volume of 0.27 Mton has potential 
of resolving neutrino mass hierarchy up to sin^ 20i3 > 0.03 (0.055) at 2cf (Sct) CL for any values 
of 6 and at the same time has the sensitivity to CP violation by 4 + 4 years running of and i>e 
appearance measurement. The significantly enhanced sensitivity is due to clean detection of mod- 
ulation of neutrino energy spectrum, which is enabled by cancellation of systematic uncertainties 
between two identical detectors which receive the neutrino beam with the same energy spectrum 
in the absence of oscillations. 
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I. INTRODUCTION 



Despite the great progresses made in varying neutrino oscillation experiments 0, S 
0, S| which revealed the structure of the lepton flavor mixing matrix, the Maki-Nakagawa- 
Sakata (MNS) matrix there still remains several unanswered questions. One of the most 
intriguing one is the problem of neutrino mass hierarchy, a tantalizing dualism of the normal 
ordering and the inverted one where the former (latter) implies that a pair of neutrinos 
with a mass gap responsible for solar neutrino oscillations are lighter (heavier) than the 
third neutrino. Picking out one from the two alternatives of neutrino mass patterns may 
bring us one of the most significant hints for underlying physics of neutrino mass. Another 
challenging goal of the next-generation neutrino oscillation experiments is to uncover leptonic 
CP violation. It will not only shed light to deep structure behind the complete parallelism 
between quarks and leptons, but also may give us a hint for understanding baryon number 
asymmetry in the universe [7|. 

In this paper, we develop a new experimental strategy of resolving the above two un- 
knowns at the same time. For the sake of concreteness, we confine ourselves into a general 
framework of the neutrino program supported by Japan Proton Accelerator Research Com- 
plex (J- PARC) located in Tokai village 0- In particular, we focus on its phase II project 
with upgraded beam of 4 MW and a megaton (Mton) water Cherenkov detector, Hyper- 
Kamiokande (HK). The fiducial mass of this detector, in its current design, is 0.54 Mton. 
We propose, instead of placing a 1 Mton detector in the Kamioka mine, to break it into 
two 0.5 Mton detectors^ and place the one in Kamioka and the other one at somewhere in 
Korea.^ Use of two identical detectors has a definite advantage that the systematic errors 
largely cancel between detectors, which allows clean measurement of modulation of neutrino 
energy spectrum due to oscillations. 

Principle of determination of mass hierarchy is quite simple. It may be carried out by 
measuring interference between vacuum and matter effects in neutrino oscillations. However, 
it requires sufficiently long baseline of the order of ~ 1000 km. Then, the natural possibility 
one can think of is to place a second detector alon g th e beam direction from Tokai to Kamioka 
at somewhere in Eurasian Continent 11, ^ 3, 3- Despite many attempts in developing 



ideas for determining mass hierarchy along the line quoted above, no immediately feasible 
proposal has emerged so far. Probably, the closest one to the goal is the NOz/A project 
(with an upgraded proton driver) in the United States jl5l |. 

It is expected in phase II of the J- PARC neutrino program that the leptonic CP violating 
phase 6 can be measured, if ^13 is not too small, by combining 2 and 6 year measurements of 
appearance channels z/^ Vc and z/^ — > z/g, respectively 0. Though powerful with enormous 
statistics enabled by intense neutrino beam and a huge detector, the measurement of the 
leptonic CP violating phase 5 in the T2K experiment suffers from the obstruction due to the 
unknown mass hierarchy, the fact observed with varying degree of robustness and contexts 
3 3 XL) 3l- ■'■^ ^ nutshell, CP violating solution can be confused with CP conserving one 
due to the degeneracy associated with the sign of ^m\^ l2], i.e., the undetermined mass 
hierarchy. Thus, the problems of measuring 5 and determining mass hierarchy inherently 



^ A two-detector method for measuring CP violation is first described in 
^ The present proposal may not require significant modification of the design of HK, because the current 

design of HK consists of two tanks (detectors) located in two separated caves due to the constraint on the 

size of the available site [lo| . 
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couple with each other through the parameter degeneracy. 

In this paper we explore a completely different strategy from these previous attempts. 
The tightly coupled feature of problems of determining mass hierarchy and 6 naturally 
suggests a new radical approach of simultaneous determination of these two quantities. For 
this purpose, we focus on a two detector system, one detector in Kamioka and another one 
at somewhere in Korean peninsula. 

The T2K experiment will use an off-axis beam, where the direction of the charged pions 
(with some contamination of kaons and muons) is 2 to 3 degrees different from the direction 
to the neutrino detector from the target. By selecting an angle between the directions of 
the charged particles and the neutrinos (off-axis angle) low energy neutrino beam whose 
maximum energy ranges from about 500 MeV (for 3 degree off-axis angle) to 800 MeV (for 2 
degree off-axis angle) can be produced. Because of the curvature of the surface of the earth, 
the direction to Kamioka, which is 295 km away from the target, is 1.3 degrees below the 
horizon as seen from the neutrino production target. According to the design of the T2K 
beam-line, the 2-to-3-degree off-axis beam is realized by directing the charged pion beam 
by 1.0 degree north of HK and 3.0 to 4.1 degrees below the horizon. The Korean peninsula 
is about 1000 to 1250 km away from the J-PARC facility. Therefore, the direction to the 
Korean peninsula is 4.5 to 5.6 degrees below the horizon. With this beam-line configuration, 
we find that the off-axis beam with the off-axis angle larger than 1 degree is automatically 
available in South Korea 0| for the 2.5 degree off-axis beam in Kamioka. We, therefore, 
find that the 2.5 degree (or more generally any angle between 2 and 3 degrees) off-axis 
beam in Kamioka and Korea is simultaneously available with the present T2K beam-line 
configuration. Later we find that the identical beam spectrum at two different positions is 
very powerful in determining the mass hierarchy. 

In Sec. m we develop our strategy of how to resolve the mass hierarchy and CP degen- 
eracies. In Sec. mil we present a pedagogical discussion of the parameter degeneracy to a 
minimal level for reader's understanding, and then illustrate how it can be resolved by spec- 
tral informations. Those who are familiar to the phenomenon should skip this section. In 
Sec. IIVI we present quantitative analysis of sensitivity of resolving the parameter degeneracy 
and determining neutrino mass hierarchy. In Sec. |3 we give concluding remarks. 



II. RIGHT AND WRONG WAYS OF USING THE KOREAN DETECTOR 



A. A detector in Korea is not always powerful 



We start by giving a cautionary remark. One must be aware that a detector in Korea 
which receives the same neutrino beam from J-PARC facility is not immediately a good idea 
for detecting the matter effect and thereby determining the neutrino mass hierarchy. To 
understand the point, we draw in Fig. ^regions spanned by probabilities of neutrino and 
anti-neutrino appearance probabilities Pii^fi —>■ i^e) and -P(i^^ —>■ i^e) in a bi-probability space 
"il. For a given ^13 variation of 5 from to 27r gives an ellipse and it forms a region when 



^13 is varied. Left and right plots in Fig. ^ are, respectively, for a detector in Kamioka and 
a one m Korea. (See for a detailed explanation of properties of the bi-probability plot.) 

The two regions with dark-gray (red) and light-gray (blue) correspond to the normal 
(Amg^ > 0) and inverted (Amg^ < 0) mass hierarchies, respectively. Here, we define the 
mass squared difference of neutrinos as Am|^ = m'j — mf where rrii is the eigenvalue of the 
ith mass-eigenstate. As one can see, the fraction of overlap between the two regions in the 
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Korean detector is as large as the one in Kamioka, indicating that there exists confusion of 
mass hierarchy in Korea as heavy as that in Kamioka. If one rehes on the usual method of 
comparing P(z/^ — >■ z/e) to Pii'/i — >■ I'e), the Korean detector is as bad as that in Kamioka. 
Of course, it implies that overall the Korean detector is worse because the number of events 
is an order of magnitude smaller for a given flux. 

Kamioka (L=295 km) Korea (L= 1 050 km) 




P(v^->v^) [%] 

FIG. 1: Regions spanned by neutrino and anti-neutrino appearance probabilities /'(f^ fc) and 
P(z/^ in a bi-probability space when CP phase 6 and ^13 are varied. Neutrino energy is 

taken as 650 MeV, and the baseline distances to Kamioka and Korea detectors are 295 km and 
1050 km, respectively. The two regions with dark-gray (red) and light-gray (blue) correspond to 
the normal {Am^i > 0) and inverted (Am^^ < 0) mass hierarchies, respectively. The left and right 
panels are for a detector in Kamioka and in Korea, respectively. A single filled circle in the left 
panel corresponds to four degenerate solutions allowed by the measurement in Kamioka, and they 
split into four different points in the right panel, which may be resolved by the Korean detector. 
Throughout this paper, mixing parameters other than ^13 and 6 are fixed to be Am|]^ = 7.9 x 10~^ 
eV^, sin^6'i2 = 0.29, Am|]^ = ±2.5 x 10^^ eV^ and sin^ 26*23 = 1, unless otherwise stated. Matter 
densities along the neutrino trajectories are assumed to be constant as 2.3 and 2.8 g/cm^ for 
Kamioka and Korea, respectively. 

The behavior we just saw above can be understood qualitatively by looking into the 
appearance oscillation probability P{i^^ i^e) (and its charge-conjugation) valid to first- 
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order in matter effect [l9l |. 
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In a(x) = y^GpNeix) [20] where Gi;' is the Fermi constant, Ne{x) denotes the electron 
number density at x in the earth, J^. (= 012^120^3513023523) denotes the reduced Jarlskog 
factor, E is the neutrino energy, L is the basehne distance, and the upper and lower sign 
± refer to the neutrino and anti-neutrino channels, respectively. We have added a term 
due to solar mass scale oscillations, the last term in (^, whose value of order 0.01 at a 
Korean detector is non- negligible. One can easily observe that when we move from the 
first to the second oscillation maxima, the size of CP phase effects become larger by a 
factor of 3 while the matter effect stays the same. At off the oscillation maxima the first 
term in last line of makes additional contributions which render the fraction of matter 
effect larger, but only up to a level given in Fig. ^ due to a modest size of the matter effect, 
aL = 0.54(p/2.8 g cm~^)(L/1000 km). Therefore, a second detector at the second oscillation 
maximum, though attractive because of a factor of 3 larger effect of CP phase 6 is not the 
best place to discriminate the mass hierarchy [l3|. 

In summary, it appears that there is no obvious merit of placing a detector in Korea to 
measure the matter effect. To really utilize the attractive feature of a factor of 3 larger CP 
effect we must go beyond the level of our above discussions, as we will pursue in the next 
subsections. 



B. Energy dependence is far more dynamic if seen by the Korean detector 

To uncover possible advantages of the Korean detector we have examined how the appear- 
ance probabilities depend upon the energy and baseline. We observed that after neutrinos 
pass through the second oscillation maximum the appearance probabilities sharply fall. The 
behavior, together with an enhanced matter effect, produces a dynamic behavior of the en- 
ergy spectrum of oscillated neutrinos. The features of the oscillation probabilities suggest 
the baselines between the second oscillation maximum and the subsequent minimum, be- 
tween about 900 and 1200 km, as appropriate ones. It is an accidental coincidence that the 
mountainous area with which neutrino beam from J- PARC facility first encounter is located 
at about 1000-1100 km away from Tokai village. Therefore, we take the baseline distance 
of 1050 km as a typical distance to the Korean detector throughout our analysis in this pa- 
per. Noting that Aml-^L/2E = 3.1l7r(Am2^/2.5 x lO'^ eV')(L/1000 km)(E/0.65 GeV)-\ 
L = 1050 km corresponds to about 10% off the second oscillation maximum at peak energy. 
The choice matches with our earlier observation that it is better to stay somewhat away 
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from the second oscillation maximum to enhance the matter effect so that the last term in 
Eq. makes some contributions to the probabilities. 

To demonstrate the dynamical behavior of the oscillated neutrino spectrum, we draw 
another bi-probability plot in Fig. |21 to represent how different are the energy dependence 
of the oscillation probabilities when they are seen by detectors in Kamioka and in Korea. 
We notice the striking difference between the left and right panels of Fig. |21 In contrast to 
the quiet behavior in Kamioka, the ellipses are larger in size and are extremely dynamic in 
movement in Korea, when the energy is varied from 0.5 to 0.8 GeV. 

Let us summarize the notable differences between Kamioka and Korea, because they are 
all crucially important as the building blocks of our strategy. 

• The variation of the appearance probabilities due to 6 is larger by a factor of 2-3 in 
Korea compared to that in Kamioka. It is notable that the ellipses with energies from 
0.5 to 0.8 GeV span almost entire area of triangular region of appearance probabilities 
less than 0.07. It is partly due to about a factor 3 larger CP effect, and to the enhanced 
solar term which scales approximately as (see Eq. (jT])). 

• The energy dependence of the ellipse is very dramatic in the Korean detector. When 
the energy is varied from 0.5 to 0.8 GeV, the slope of the major axis moves coun- 
terclockwise (clockwise) by more than 71/2 in the case of normal (inverted) hierarchy, 
which is represented with the symbol exhibited in Fig. |21 (For electronic version, the 
ellipses of the normal and the inverted hierarchies are represented by warm and cold 
colored lines, respectively.) It implies that CP phase 6 affects differently onto different 
part of the energy spectrum. Notice that positive slope of the major axis has never 
been observed in region around the first oscillation maximum. In contrast to such 
violent behavior the energy dependence of the ellipse is very quiet in Kamioka. 

This strong energy dependence of the appearance oscillation probabilities in Korea, which 
are different for normal and inverted hierarchies and are helped by their enhanced size, is 
the key to the enormous sensitivity of the two-detector complex as we will fully explore in 
Sec. US 

We want to mention here about the BNL proposal in which a wide-band beam from 
upgraded AGS will be detected by the far detector located at 2950 km to observe the effects 
of multiple (1-3) oscillation maxima jij]. The dramatic difference in the energy dependence 
of the neutrino oscillation probability between the first and the second oscillation maxima 
may have been contributed to the high sensitivity (despite such a long baseline) obtained in 
the proposal. 

C. Determination of mass hierarchy through resolving degeneracy 

Now we have reached to the point to establish our new strategy for the determination of 
the neutrino mass hierarchy. In a nutshell we do it by resolving the parameter degeneracy 
associated with the unknown sign of Amg^. While we postpone an organized discussion 
of the parameter degeneracy to Sec. IllH it will be shown there that to resolve the Arrig^- 
sign degeneracy, the spectrum shape analysis of the appearance events in two detectors in 
different baseline (or at different energy) is crucial. Then, by placing the two detectors in 
Kamioka and in Korea, we can fully utilize the remarkable difference in energy dependence 
of the oscillation probabilities to solve the parameter degeneracy. 
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Kamioka (L=295 km) 



Korea (L=1050 km) 




P(v -v; [%] 



FIG. 2: Energy dependences of the oscillation probabilities for sin^ 2^13 = 0.05 are represented 
by plotting ellipses (which results as 6 is varied from to 2tt) in bi-probability space for various 
neutrino energies from 0.5 to 0.8 GeV. Other mixing parameters are fixed as in the caption of 
FigH The left and the right panels are for detectors in Kamioka and in Korea, respectively. The 
correspondence between energies and the line symbols are denoted in the figure; The ellipses in 
upper 4 symbols (warm colors) indicate the ones of normal mass hierarchy (Ani^^ > 0) and the 
one of lower 4 symbols (cold colors) the ones of inverted mass hierarchy (Am^^ < 0). 



There exists great merit of using two identical detectors together with the neutrino beam 
of identical spectrum shape (in the absence of oscillations) which makes demonstration of 
the CP violation and resolution of the mass hierarchy degeneracy much cleaner.^ We will 
also fully utilize this characteristics of Kamioka-Korea two identical detector complex in our 
analysis in Sec. IIVI 



III. PARAMETER DEGENERACY AND ITS RESOLUTION BY SPECTRAL 
ANALYSIS 

In this section, we briefly summarize the parameter degeneracy in neutrino oscillations 
[1^ H^, 113, Ho, H^l to the extent that is necessary to understand the rest of this paper. 
We then illustrate how it can be resolved with use of the Kamioka-Korea two identical 
detectors, which greatly differ in energy dependences they observe. Those who are familiar 



^ We are not aware to whom the idea of error cancellation with use of identical detectors can be attributed. 
It was utilized, for example, in the Bugey experiment ji^j and there were many examples of constructing 
near detectors more or less identical with far detectors. 



7 



to the parameter degeneracy should skip this section and directly go to Sec. IIVI where full 
details of the quantitative analysis of resolving power of the degeneracy will be presented. 

It should be obvious from Fig.[T]that a set of measurement of -P(t'^ — > t'e) and P{i^n I'e) 
at a particular baseline and an energy allows four solutions (crossing points of the eUipses) of 
^13 and 6. The origin of the four-fold solutions may be understood as the intrinsic degeneracy 
of 6^13 and 5 (a single crossing point allows two ellipses), duplicated by the unknown sign 
of Amg^ Jl^ which admits opposite Amg^-sign solutions. The mixed Amg^-sign degenerate 
solution occurs in the cross-hatched region in Fig. ^ A degenerate solution in Kamioka (left 
panel) splits into four separate points (right panel) in the Korean detector. The importance 
of the Korean detector is clearly understood. 

In fact, the four-fold degeneracy is further multiplied by the octant ambiguity of ^23 for 
a given sin 26*23 iS^Jj leading to the total eight-fold degeneracy But we do not try to 
resolve the ^23 degeneracy in this paper by simply assuming that 6*23 = 7r/4 in most part of 
our analysis. We, however, examine the effect of non-maximal ^23 on resolution of neutrino 
mass hierarchy and CP degeneracy.^ 

Now we present in Fig. |21a four- fold degenerate solutions to show how it looks like in the 
plane spanned by sin^ 2^13 and 5 obtained as a result of our analysis whose details will be 
fully explained in Sec. IIVI These solutions are obtained by assuming the operation of the J- 
PARC experiment for total 8 years running with and z/e appearance channels. The upper 
figure: 2 years running of the neutrino mode and 6 years running of anti-neutrino mode with 
single detector of 0.54 Mton of the fiducial volume placed in Kamioka only. The lower figure: 
4 years of running of both in neutrino and anti-neutrino modes with two detectors of 0.27 
Mton of the fiducial volume placed one in Kamioka and the other in Korea. In each figures 
the left and the right panels are for cases with true values of ^13 given by sin^ 2^13 = 0.02 
and sin^ 26*13 = 0.005, respectively, with positive value of Amg^, as inputs. 

Let us understand the features of the four solutions represented in Fig. El First, we focus 
on the right panel of Fig. 01 We observe two distinct pairs of solutions, the same Amg^^-sign 
(intrinsic) degeneracy and the mixed Aml^-sign degeneracy. 

• the same Am|i-sign degeneracy 

One can identify a pair of the solutions in each mass hierarchy, which obviously corre- 
sponds to the intrinsic degeneracy. It is known that the two solutions of 5 are approximately 
symmetric with respect to reflection at 5 = 7r/2, that is ^2 — (tt — 5i) < 0.027r [27,,^]. Within 
the same mass hierarchy, the transformation 5 ^ (vr — 5) produces flipping the sign of cos 5 
term in the appearance probability in Eq.(0). Therefore, the oscillation probabilities with 
these two degenerate solutions have different energy dependence, and hence they are discrim- 
inable by doing spectrum analysis. This is what happens, though incompletely, in the case 
of sin^ 2^13 = 0.005 (right panels). For sin^ 2^13 = 0.02 (left panels) the spectral information 
is powerful enough to resolve (almost) the intrinsic degeneracy even for a single detector 
in Kamioka, as seen in the left panels of the upper plots in Fig. |21 Use of the spectral 



^ It is well recognized that it is very difficult to resolve the degeneracy associated with ^23 if we rely only 
on accelerator measurement using beam. It in turn implies that our discussion in this paper may be 
stable against inclusion of deviation of ^23 from the maximal. Lifting the ^23 degeneracy may be done 
only when the accelerator measurement of v^^ — )■ is combined with cither pure measurement of ^13 by 
reactor experiments or complementary information from the silver channel, Vr |2^. 
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informations to solve the parameter degeneracies has been discussed by many authors; See 



solve tne_p£ 



for example |12l . [1 

• the mixed Amg^-sign degeneracy 

One can also identify another pair of solutions which are still allowed even at 68% CL 
in the opposite-sign Am^i panels. This is the degenerate solution due to the unknown sign 
of Amg]^ (or type of mass hierarchy). It arises because of the approximate symmetry under 
the simultaneous transformation (which applies if the matter effect is small) 

6 ^ IT — 6 (mod.27r), 

Ami, ^ -A^3i- (2) 

Since it is valid apart from tiny corrections of the order of sin^ 29i3{Am2i/ Amli) these two 
solutions cannot be distinguished only by spectrum analysis at the single detector with small 
matter effect. It is the reason why we need a Korean detector to resolve the degeneracy due 
to mass hierarchy. The fact that the same transformation of 6 appears in Q explains the 
reason why two pairs of the mixed Amg^-sign solutions forms X (cross) shape in Fig. |21 

To have a feeling on how the spectrum information helps to resolve the degeneracies, we 
present in Fig. |3 the energy dependences of the oscillation probabilities to be observed at 
detectors in Kamioka (left panels) and in Korea (middle and right panels). In the right 
panels, we artificially switch off the matter effect to illuminate its role in resolving the 
degeneracy. We pick up a set of mixing parameters which correspond to a degenerate 
solution in Kamioka from the rate only analysis, which are shown in the figure, and indicate 
how these four solutions differ in energy dependences in Kamioka and Korea. The value 
taken for ^13, sin^ 2^13 ~ 0.02, is large enough for the Kamioka-Korea two-detector complex 
to solve all the degeneracies, but is too small for a single detector in Kamioka to do the job, 
as indicated in the left panels of Fig. |21 It is obvious from Fig. |3] that the Korean detector 
sees very different energy spectrum of neutrinos. We also note the several characteristic 
features: 

(1) In Kamioka (left panels), there are some differences between the same-sign degeneracy 
and the mixed-sign degeneracy. In the former case energy dependences of a pair of the 
solutions are somewhat different, leading to resolution of intrinsic degeneracy (as is seen in 
the left panels of Fig. ^ if the statistics is sufficient and the systematic error is controlled. 

(2) In Korea (middle panels), the energy dependences of the four solutions are so distinct, 
and it is conceivable that the spectral information solves all the four-fold degeneracies, 
provided that the statistics is sufficiently large. 

(3) Comparison between the middle and the right panels of Fig. |3 reveals that the matter 
effect plays a key role in resolving the Aml^-sign degeneracy (see solid red and blue dashed 
curves) . 



IV. SENSITIVITIES OF RESOLVING MASS HIERARCHY AND CP DEGENER- 
ACY 

In the first subsection we give a full details of our analysis, including treatment of ex- 
perimental errors, treatment of background, and the statistical procedure which is used to 
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FIG. 3: Examples of the parameter degeneracy. The true solutions are assumed to be located at 
(sin^ 2^13 and 6) = (0.005, 7r/4) (left panels) and (0.02, 7r/4) (right panels) with positive sign of 
Am^^, as indicated as (green) stars. Three contours in each figure correspond to the 68% (dotted 
lines, blue), 90% (solid lines, black) and 99% (dashed lines, red) C.L. sensitivities, which are defined 
as the difference of the being 2.30, 4.61 and 9.21, respectively. The upper and lower panels are 
for a detector in Kamioka and for two detectors one in Kamioka and the other in Korea (twin HK), 
respectively. 
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FIG. 4: Neutrino oscillation probabilities corresponding to a four-fold degenerate solutions obtained 
by measurement in Kamioka by the rate only analysis are plotted as a function of neutrino energy. 
Left panels: appearance probabilities in Kamioka. Middle panels: appearance probabilities in 
Korea. Right panels: appearance probabilities in Korea, but without the matter effect. 



investigate the sensitivity of resolving the degeneracies. In the next subsections we give our 
results and perform stability checks of the results against the uncertainty of the matter den- 
sity, the change in the systematic errors, and the possible deviation of from the maximal 
value. 



A. Assumptions and the definition of 

In order to understand the advantage of the two detector system at 295 km (Kamioka) 
and 1050 km (Korea), we carry out a detailed analysis. We use various numbers and 
distributions available from references related to T2K [s^ • Here we summarize the main 
assumptions and the methods used in the analysis. 

We use the reconstructed neutrino energy for single- Cherenkov-ring electron-like events 
(about 80% of which are quasi-elastic charged-current (CC) Pe interactions). The resolution 
in the reconstructed neutrino energy is 80 MeV for quasi-elastic events. We assume that 
sin^ 26'23 and | Am|^ | should be known precisely by the time when the experiments we consider 
in this report will be carried out. We take for these parameters sin^2^23 = 1 and Amg^ = 
±2.5xlO-3eV^ Hence, we assume that the energy spectrum of the beam is the one expected 
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by the 2.5 degree off-axis-beam in T2K. The shape of the energy spectrum for the anti- 
neutrino beam is assumed to be identical to that of the neutrino beam. The event rate for 
the anti- neutrino beam in the absence of neutrino oscillations is smaller by a factor of 3.4 
due mostly to the lower neutrino interaction cross sections and partly to the slightly lower 
neutrino flux. 

In the T2K experiment, 28 background events are expected for the reconstructed neutrino 
energies between 350 and 850 MeV for (0.75 x 0.0225 x5)MW-Mton-yr measurement with the 
neutrino beam. The energy dependence of the background rate and the rate itself are taken 
from [32] . The background rate is expected to be higher in the lower neutrino energies. The 
expected number of signal events is 122 for sin^ 2^13 =0.1 with the same detector exposure 
and beam, assuming the normal mass hierarchy and 6 = 0. The background rate for the 
anti-neutrino beam is smaller than that for the neutrino beam by a factor of 1.6. Since the 
signal rate for the anti-neutrino beam is smaller by a factor of 3.4, the signal to noise ratio 
is worse for the anti-neutrino beam than that for the neutrino beam by a factor of about 2. 

We assume that the experiment has a front detector that measures the background and 
the signal detection efficiency. We assume that both the overall background rate and the 
background shape (i.e., the energy dependence of the background) are understood within 
the uncertainty of 5%. More details on the assumed uncertainty in the background will be 
discussed later in this subsection when we define the for the analysis. We stress that in 
the present setting the detectors located in Kamioka and in Korea are not only identical but 
also they receive neutrino beams with the same energy distribution (due to the same off-axis 
angle) in the absence of oscillations. Therefore, the dependence on the assumed value of 
the experimental systematic errors must be weak, as we will verify by repeating the same 
analysis with different values of the uncertainties. As we will see later, the error cancellation 
between identical detectors gives us a very powerful constraint in our analysis. 

We compute neutrino oscillation probabilities by numerically integrating neutrino evolu- 
tion equation under the constant density approximation. The average density is assumed 
to be 2.3 and 2.8 g/cm^ for the matter along the beam line between the production target 
and Kamioka and between the target and Korea, respectively.^ We assume that the number 
of electron with respect to that of nucleons to be 0.5 to convert the matter density to the 
electron number density. 

Fig. El shows an example of the energy spectrum of electron events to be observed in 
Kamioka and Korea for 4 years of neutrino beam and 4 years of anti-neutrino beam. Even 
with the 4 MW proton beam and with 0.27 Mton fiducial volume, the number of expected 
events in Korea is limited due to the large distance between the production point and the 
detector. However, it is seen that the overall signal rate and the energy dependence of the 
signal events are different between Kamioka and Korea. Therefore, we use these features to 
obtain information on neutrino oscillation parameters. 

The statistical significance of the measurement considered in this paper was estimated 



^ The former number, as quoted in [33, is the average value of the matter density profile which is obtained 
by M. Komazawa at Geographical Survey Institute in Japan by measuring surface gravity in many points 
around the path from KEK to Kamioka for the K2K experiment. The latter one is a commonly taken 
value for the matter density in earth crust, which is roughly equal to the average density of the distribution 
expected by "Preliminary Reference Earth Model" [sj]. 
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FIG. 5: Examples of electron events to be observed in Kamioka and in Korea for 4 years of 
operation of the neutrino and anti-neutrino beams are presented as a function of reconstructed 
neutrino energy. The fiducial masses are taken to be 0.27 Mton for both the detectors in Kamioka 
and in Korea. The open circles show the background events. The solid circles and dashed lines show 
the expected energy spectrum of signal events with and without the detector effects, respectively, 
sin^ 2^13 =0.1, 6 = 7r/2 and normal mass hierarchy are assumed in simulating the events. 



by using the following definition of j35| : 

x-E(E«^«^).E(f)^ (3) 

k=l \i=l J j=l ^^i^ 

N{err = A^f ^ ■ (1 + E ^^■) + ^^"^^ ■ (1 + /3 ■ ^3) . (4) 

The first sums are for the number of observed single-ring electron events in the i^^ energy 
bin, A^(e)°'^^ is the number of events to be observed for the given oscillation parameter set, 
and N[eY^'^ is the expected number of events for the assumed sin^26'i3, 5 and the mass 
hierarchy in the analysis, where k = 1,2,3 and 4 correspond to the four combinations 
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of the detectors in Kamioka and in Korea with the neutrino and anti-neutrino beams, re- 
spectively. Both N{e)°^^ and A^(e)^^^ include background events. The energy ranges of the 
five energy bins are respectively 400-500 MeV, 500-600 MeV, 600-700 MeV, 700-800 MeV 
and 800-1200 MeV. cTj denotes the statistical uncertainties in the expected data. The energy 
resolution of 80 MeV is taken into account with a Gaussian resolution function. 

During the fit, the values of N^^^ are recalculated to account for neutrino oscillations, 
and so are the systematic variations in the predicted rates due to the uncertainties in the 
estimated signal and background. The overall background normalization is assumed to be 
uncertain by ±5%, and the effect is taken into account through ei in Q. In addition, it is 
also assumed that the background has an energy dependent uncertainty with the functional 
form of {{E^{rec) — 800 MeV)/400 MeV) x (1 + €2). The energy- dependent part is also 
assumed to be uncertain by 5%. The uncertainty in the detection efficiency of the electron 
and positron signals is assumed to be 5%. In summary, dj = 0.05 for j = 1,2, and 3. 

N^"-^ is the number of background events for the i^^ bin. We note that the number 
of the background events at detectors in Kamioka and in Korea are related simply by 
{L Korea / LKamiokaY 1 where L is the distance between the neutrino production point and 
the detector, f^^'^s^ai ^-^^ number of events that are appeared by neutrino oscillations, 
and depend on sin^2^i3, 5 and the mass hierarchy. The uncertainties in Nj^'^ and ]\fp9nai 
are represented by 3 parameters ej. The parameter /j represents the fractional change in 
the predicted event rate in the i*^ bin due to a variation of the parameter ej. The third 
sum in the definition collects the contributions from variables which parameterize the 
systematic uncertainties in the expected number of background events. During the fit, these 
3 parameters are varied to minimize for each choice of the oscillation parameters. 

B. Sensitivity achievable by two-detector complex 

We have shown in Fig. El typical results of the sensitivity analysis by taking as an example 
two values of 613, sin^ 26*13 = 0.005 and 0.02 with normal mass hierarchy as inputs. In the 
upper and the lower figures the cases of a single 0.54 Mton detector in Kamioka and two 
detectors of 0.27 Mton each in Kamioka and in Korea, respectively, are examined. Now, 
we proceed to the sensitivity analysis of mass hierarchy and CP violation, which covers the 
entire region of the relevant parameter space. 

We examine five different detector configurations, namely 5 cases of volume ratios of 
Kamioka to Korean detectors, keeping the total volume constant; (a) a single 0.54 Mton 
fiducial mass detector (the total mass of which should be about 1 Mton) in Kamioka, (b) 
0.27 Mton detectors in Kamioka and in Korea, (c) a 0.16 Mton detector in Kamioka and a 
0.38 Mton detector in Korea (volume ratio of 3:7), (d) 0.05 Mton detector in Kamioka and 
a 0.49 Mton detector in Korea (volume ratio of 1:9), and (e) a single 0.54 Mton detector in 
Korea. It is assumed that the experiment will continue for 8 years with 4 year of neutrino 
beam and 4 years of anti-neutrino beam. We use the same definition of as in Q which 
embodies cancellation of errors between identical detectors even in the cases of two detectors 
with different volumes. For this simplification and the other reasons our analysis cannot be 
considered as a real optimization procedure for the volume ratio, but we feel that such study 
may give some insights on how the sensitivity varies with the detector mass ratio. 

We show in Fig. IHl the contours of sensitivity to the mass hierarchy at 2 and 3 stan- 
dard deviations on 5-sin^ 26*13 plane for these 5 detector configurations. Here, the 2 and 3 
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FIG. 6: 2(thin lines) and 3(thick lines) standard deviation sensitivities to the mass hierarchy for 
the volume ratio of Kamioka to Korea 1:0 (0.54 Mton fiducial mass detector in Kamioka, indicated 
by dashed lines, blue), the volume ratio 1:1 (0.27 Mton detectors in Kamioka and in Korea, solid 
lines, black), the volume ratio 3:7 (a 0.16 Mton detector in Kamioka and a 0.38 Mton detector 
in Korea, dash-dot lines, yellow), the volume ratio 1:9 (a 0.05 Mton detector in Kamioka and a 
0.49 Mton detector in Korea, dash-dot-dot lines, green), and the volume ratio 0:1 (a 0.54 Mton 
detector at Korea, dotted lines, red). 4 years running with neutrino beam and another 4 years 
with anti-neutrino beam are assumed. 



standard deviations are defined to be |Xmm("^^o^g hierarchy) — x^j„(true hierarchy)! > 4 
and 9, respectively. In the upper and the lower panels of Fig. IHl cases of the normal and 
the inverted mass hierarchies as a true mass hierarchies are presented. As expected, better 
sensitivities are obtained for detector configurations with detectors both in Kamioka and in 
Korea. While the best sensitivity is obtained for the case of (c) 3:7 volume ratio, we find 
that the sensitivity on the mass hierarchy does not depend strongly on the mass ratio of the 
two detectors. It is interesting to note that there are some sensitivities to mass hierarchy 
in a limited region of 6 even in the setting with a single detector in Kamioka, as expected 
110. 

We show in Fig. [7| the sensitivity to leptonic CP violation by drawing the contours above 
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FIG. 7: The same as in Fig. IHlbut for the sensitivity to leptonic CP violation. The correspondence 
between line symbols and detector configurations is the same as in Fig. |H1 

which hypothesis of 5 = or vr can be rejected at 2 and 3 standard deviations. The 2 and 
3 standard deviations are defined to be IXmini^ 7^ or vr) — x^j„(true value of 6)\ > 4 and 
9, respectively. As in the upper and the lower panels of Fig. the cases of the normal and 
the inverted mass hierarchies as a true mass hierarchies are presented. Notice that while 
the type of true hierarchy is assumed in each case, the fit is performed without assuming 
the mass hierarchy. This statement applies to all the figures presented in this section. For 
small sin^ 2^13 a single 0.54 Mton detector in Kamioka has the best sensitivity. Whereas for 
large sin^ 26*13, the two detector configuration gives the better sensitivity. It is due to the 
fact that in the two detector configurations, the mass hierarchy can be resolved for large 
sin^ 26*13, and therefore the value of 6 is uniquely measured in the case of the two detector 
configurations. Among the two detector configurations, 0.27 Mton detectors at Kamioka 
and Korea gives the best sensitivity to the CP violation measurement. 

We observe that the sensitivities to the mass hierarchy and CP violation do not change 
appreciably as the volume ratio is varied. A larger Korean detector seems to be favored 
to resolve the mass hierarchy. On the other hand, the sensitivity to CP violation disfavors 
the larger Korean detector, and in fact the best choice would be a 0.54 Mton detector in 
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Kamioka for small sin^2^i3 values. Overall, two 0.27 Mt detectors one in Kamioka and the 
other in Korea seems to be close to the optimal choice. Anyhow the difference is so small 
that the real optimization process must involve other factors, such as the site availability 
and conditions, as well as additional physics capabilities, etc. 

We note, in passing, that in the case of a single 0.54 Mton detector in Kamioka the spectral 
analysis is so powerful that allows larger sensitivity region of discovering CP violation than 
the one presented in Q (which relied on rate only analysis), as one can see in Fig. [7| It is 
highly nontrivial because we used relaxed values of systematic errors of 5%. Moreover, one 
can recognize the power of the spectrum information to solve the intrinsic degeneracy. 

iFrom the results obtained in this subsection, we select out the two identical detector 
configuration with 0.27 Mton each in Kamioka and in Korea for further studies of the 
sensitivities of resolving the mass hierarchy and the CP degeneracy, which is to be carried 
out in the following subsections. 

C. Optimization of neutrino and anti- neutrino running time 

Since we found that the best detector configuration is the identical detectors in Kamioka 
and in Korea, we search for the best option for dividing assumed whole period of 8 years into 
the neutrino and anti-neutrino running. We examine the following three cases: (i) 8 years 
of the neutrino beam, (ii) 4 years of the neutrino and anti-neutrino beams (4+4), and (iii) 

2 years of the neutrino beam and 6 years of the anti- neutrino beam (2+6). The last option 
is studied in Ref. where the sensitivity on the CP violation is discussed with a single 
detector in Kamioka. Figure |H1 (upper panels) shows the sensitivity to the mass hierarchy in 
5-sin^ 2^13 space at 2 and 3 standard deviations for the three beam options. We find that the 
(4+4) option gives the best sensitivity while the (2+6) option gives similar but somewhat 
lower sensitivity. 

Figure |H1 (lower panels) shows the sensitivity to the leptonic CP violation in 5-sin^ 26'i3 
space at 2 and 3 standard deviations for the three beam options. The (4+4) and (2+6) op- 
tions have almost the same sensitivity to the leptonic CP violation. We note that this study 
suggests that the leptonic CP violation can be demonstrated by running the experiment 
with the neutrino beam alone for some values of sin^ 2^13 and 6 though the sensitivity is 
substantially lower than the neutrino plus anti-neutrino beam options. This can be under- 
stood because the energy spectrum data from the detectors in Kamioka and Korea contain 
information on 6. 

In Fig. El we present sensitivity regions for the mass hierarchy and CP violation at 2 and 

3 standard deviations in space spanned by sin^ 26'i3 and the fraction of 6. The plot probably 
needs explanation: For a given value of sin^ 26'i3, one can draw the region, as we did in Fig. [HI 
Fig. [71 and Fig. |H1 in the previous subsections, where the mass hierarchy can be resolved (or 
CP violation can be signaled) at 2 or 3 standard deviations. The plot can be converted to 
the one in Fig. El by computing the fraction of the region of S above the resolution contour 
at the given value of sin^ 2^13. In this way, the plot gives an alternative way of representing 
how sensitive is the experiment to the mass hierarchy or CP violation, and is used e.g., in 
the proposal of the NOi/A experiment . 

In the proposed experiment, the mass hierarchy can be resolved if sin^ 2^13 is larger than 
0.055 (0.03) at the 3 (2) standard deviation level for any 6 values. (Note, however, that the 
sensitivity is better in most of the region of 5 as can be seen in Fig. IHl) 
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FIG. 8: 2(thin lines) and 3(thick lines) standard deviation sensitivities to; the mass hierarchy 
(upper panels) and the leptonic CP violation (lower panels) for 0.27 Mton detectors in Kamioka 
and in Korea with 4 years running with neutrino beam and 4 years with anti-neutrino beam (solid 
lines, black), 8 years with neutrino beam only (dotted lines, red) and 2 years with neutrino beam 
and 6 years with anti-neutrino beam (dashed lines, blue). 
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FIG. 9: 2(thin line) and 3(thick lines) standard deviation sensitivities to mass hierarchy (left 
panel) and CP violation (right panel). Solid and dashed black lines show the cases of two identical 
detectors with fiducial mass of 0.27 Mton in Kamioka and in Korea with 4 years of neutrino and 
4 years of anti-neutrino beam running. Dotted and dash-dotted blue lines show the cases for a 
single detector with fiducial mass of 0.54 Mton in Kamioka with 2 years of neutrino and 6 years 
of anti-neutrino beams. Solid (dotted) and dashed (dash-dotted) lines show the cases for positive 
and negative mass hierarchies, respectively. 

D. Systematic uncertainties and the sensitivity 

In order to study the robustness of the results, we carry out several tests. First, we test 
stability of the results by varying the systematic uncertainties in the background estimation 
and the signal detection efficiency. We examine three values, 2, 5, and 10%. FigurelTOk shows 
the sensitivity region to the mass hierarchy in 5-sin^ 2^13 space at 2 and 3 standard deviations 
for the three values of systematic uncertainties. It is remarkable that the dependence of the 
sensitivity to the systematic error is extremely weak. We point out that, because of the 
identical energy spectrum and the detectors in Kamioka and in Korea, their background 
rates must be related simply by (i^Korea/-Z^Kamioka)^- Because of this relation, the difference 
in the signal events in Kamioka and Korea cannot be explained by the uncertainty in the 
background, and therefore the measurement of the background rate by front detectors with 
a very high precision is not crucial in the identical two detector setup. 

The matter density along the neutrino beam is not precisely known. Uncertainty in the 
matter density between J-PARC and Korea could be non-negligible, while the matter effect 
could be more important. We study the systematic effect due to the uncertainty in the 
matter density by changing only the mean matter density between the target and Korea. 
The true mean matter density is assumed to be either 2.5, 2.8 or 3.1 g/cm^, while the analysis 
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FIG. 10: (a) (left panels): 2(thin line) and 3(thick lines) standard deviation sensitivities to mass 
hierarchy for systematic uncertainties of 2% (dashed lines, blue), 5% (solid lines, black), and 10% 
(dotted lines, red) in the estimated signal and background, (b) (middle panels): the same as in 
(a) but with varying true values of mean matter density, 2.5 (dashed lines, blue), 2.8 (solid lines, 
black), and 3.1g/cm'^ (dotted lines, red). The analysis is carried out assuming the matter density 
of 2.8 g/cm^. (c) (right panels): the same as in (a) but with varying true values of sin^^23) 0.45 
(dashed lines, blue), 0.5 (solid lines, black), and 0.55 (dotted lines, red). The analysis is carried 
out assuming sin^ 623 = 0.5 

is carried out assuming the density of 2.8g/cm'^. Figure (TUb (middle panel) shows contours 
of the sensitivity to mass hierarchy in 5-sin^ 26*13 space for three cases of the matter density. 
As indicated in the figure, there is only a small difference in the sensitivity to the mass 
hierarchy. 

Even if sin^ 26*23 is measured to an accuracy of 1%, the uncertainty in sin^ 6*23 is relatively 
large. To examine the dependence of the sensitivity to 6*23, we have assumed that the 
true value of sin^ 6*23 is either 0.45, 0.5 or 0.55, while the analysis is carried out assuming 
sin^ ^23 =0.5. Figure ITUh (right panel) shows the sensitivity to the mass hierarchy in 6- 
sin^ 26*13 space for the three values of sin^ 2^23 • Although the sensitivity to mass hierarchy is 
somewhat influenced by the uncertainty of sin^ ^23, the mass hierarchy can still be resolved 
without any significant overall change in the sensitivity. The effect of the uncertainty in 
sin^ ^23 to the sensitivity to the leptonic CP violation is very small. 

In summary, we find that the conclusion we presented in the previous subsection does 
not depend strongly on the assumptions we have made. 

V. CONCLUDING REMARKS 

In this paper, we have explored physics capabilities of the two identical megaton-class 
detector complex (twin HK), one placed in the Kamioka mine and the other in Korea 
under the assumption that they receive the same neutrino beam from J-PARC with 4MW 
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beam power. Unlike the foregoing analyses of such two detector setup, we have employed 
a radical new strategy of resolving the Amg^^-sign degeneracy to determine simultaneously 
the neutrino mass hierarchy and the CP violating phase 6. We have demonstrated that the 
two-detector complex can determine neutrino mass hierarchy down to sin^ 2^13 > 0.02 (0.05) 
for ~80% coverage of the whole range of 6, and sin^ 26'i3 > 0.03 (0.054) for any value of 6 
both at 2(7 (3(t) CL, as indicated in Fig. El Importantly, the sensitivity to the CP violation 
of the current design of J-PARC phase II project is essentially kept or even enhanced at 
sin^26'i3 > 0.01. The key to the enormous sensitivities is the use of two identical detectors 
that allows significant reduction of systematic errors by cancellation. It allows unambiguous 
detection of spectral distortion of neutrino energy distribution due to oscillations, which is 
crucial to resolve the Amg^-sign (the intrinsic) degeneracy with (without) help of the matter 
effect. 

We have shown that our proposal elevates the current design of the phase II of J-PARC 
neutrino program which aims at discovering leptonic CP violation to a dual-purpose exper- 
iment, in which one can determine neutrino mass hierarchy and CP violation up to a level 
achievable by conventional superbeam experiments. We emphasize that, in our setting, the 
total volume of the two HK is the same as the planned HK with 1 Mton water in Kamioka, 
and the current design of the latter option already contains two tanks. Hence, the cost of the 
present proposal should be roughly the same as that of the conventional design. Yet, we nei- 
ther addressed the issue of site availability in Korea, nor the precise estimate of systematic 
errors which requires specification of the site as well as actual design of the detectors. 

There exist several other proposals for the next-generation long-baseline experiments in 
the world, NOz/A [l^, the very- long-baseline project j2l|, and the beta beam Js?]. With 
different assumptions and conditions, it is difficult to compare the present proposal with the 
existing ones. Nonetheless, we would like to emphasize that the present proposal is at least 
competitive to them not only on the sensitivity to the CP violation but also on resolving 
power of the neutrino mass hierarchy. 
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